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ABSTRACT 

We present deep spectroscopic observations of the classical T Tauri stars DF Tau and V4046 Sgr in 
order to better characterize two important sources of far-ultraviolet continuum emission in protoplan- 
etary disks. These new Hubble Space Telescope-Cosmic Origins Spectrograph observations reveal a 
combination of line and continuum emission from collisionally excited H2 and emission from accretion 
shocks. H2 is the dominant emission in the 1400 < A < 1650 A band spectrum of V4046 Sgr, while an 
accretion continuum contributes strongly across the far-ultraviolet spectrum of DF Tau. We compare 
the spectrum of V4046 Sgr to models of electron-impact induced H 2 emission to constrain the physical 
properties of the emitting region, after making corrections for attenuation within the disk. We find 
reasonable agreement with the broad spectral characteristics of the H 2 model, implying AT (if 2) ~ 10 18 
cm~ 2 , T{H2) = 3000^5qo° K, and a characteristic electron energy in the range of ~ 50 - 100 eV. We 
propose that self- absorption and hydrocarbons provide the dominant attenuation for H2 line photons 
originating within the disk. For both DF Tau and V4046 Sgr, we find that a linear fit to the far-UV 
data can reproduce near-UV/optical accretion spectra. We discuss outstanding issues concerning how 
these processes operate in protostellar/protoplanetary disks, including the effective temperature and 
absolute strength of the radiation field in low-mass protoplanetary environments. We find that the 
912 - 2000 A continuum in low-mass systems has an effective temperature of ~ 10 K with fluxes 
10 5 ~ 7 times the interstellar level at 1 AU. 

Subject headings: accretion — protoplanetary disks — stars: pre-main sequence — stars: individual 
(DF Tau, V4046 Sgr) 



1. INTRODUCTION 

Classical T Tauri stars (CTTSs) are characterized by 
broad Ha emission lines and ultraviolet (UV) spectra 
dominated by atomic and m olecular features t hat are 
attributed to gas-rich disks (jFurlan et al.l 120061) . The 
ages of CTTS disks (0.1 - 12 Myr; Isella et al. 2009; 
Kastner et al. 2008) indicate that the gaseous pro- 
cesses contributing to their observational characteristics 
are contemporaneous with and likely intimately con- 
nected with giant planet formation, which is thought to 
be m ostly completed on similar timescales (jAlibert et al.l 
2005). Terrestrial planet for mation is thought to occu r 
on timescales of 10 - 100 Myr (jKenvon fc Bromlevl[200fih . 
when the majority of the primordial gas is in the disk 
has been dissipated. Observations of the gas and dust 
in CTTS disks therefore probe the physical and chemical 
state of gas giant forming protoplanetary systems. 

Multi-wavelength spectroscopy is a powerful tool for 

* Based on observations made with the NASA/ESA 
Hubble Space Telescope, obtained from the data archive at the 
Space Telescope Science Institute. STScI is operated by the As- 
sociation of Universities for Research in Astronomy, Inc. under 
NASA contract NAS 5-26555. 

1 kcvin.france@colorado.edu 



making quantitative measurements of the dust and gas in 
these systems. Dust in protoplanetary disks is seen most 
clearly through the mid- an d far-infrared (IR) excess 
flux produced by warm grains (jFurlan et al.1120061) . which 
overwhelms the narrow molecular lines when observed at 
low-resolution (|Naiita et al.l I2010D . Emission from hot 
gas is produced by magnetic activity in the atmospheres 
of the centr al stars and the accretion shocks near the ste l- 
lar surface (|Johns-Krull et al.ll2000t lArdila et alJl2053 l. 
while the protoplanetary material itself can be probed 
through molecular observations of these systems. Pro- 
toplanetary disks have a multi-phase physical structure 
(see, e.g., the reviews by Woitke et al. 2009 and Dulle- 
mond & Monnier 2010) seen in molecul ar line o b serva - 
tions from the far-UV to the millimeter. iDutrevi ()2001[ ) 
reviewed the mm- wave CO measurements of outer pro- 
toplanetary disks, but understanding the complex CO 
excitation structure with radius requires additional trac- 
ers of th e bulk molecu lar material within 100 AU of 
the star (|Greavesll200l . Emission from OH, CO, C0 2 , 
H 2 0, and other biologically important species have been 
used to trace war m gas (T ~ few x 10 2 - 2000 K) 
in the inner disk (iNaiita et al.l l2003t iSalvk et al.l 120081 
IBethell fc Berginll2009ri . " 
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Metal-bearing molecules trace molecular hydrogen 
(H 2 ), the primary constituent of gas giant planets. 
However, H2 can be observed best in the far-UV 
(912 - 1650 A) bandpass, where the dipole-allowed 
molecular emission spectrum is prim arily photo-excited 
("pumped") by stel lar Lya photons (|Ardila et al.l 120021 : 
IHerczeg et alll200l . The Lya -pumping route requires 
that the second excited vibrational level (v — 2) of H2 
have an appreciable population, which requires that the 
molecules reside i n a warm (Tfflo,) > 2000 K) gas layer 
close to the star (|Herczeg et al.l [20061) . Most likely, this 
photoexcited H2 is not physically associated with the 
bulk of the colder H2 that goes into planet formation. 
The homonuclear nature of H2 means that rovibrational 
transitions are dipole forbidden, making direct detec- 
tion of the C00I-H2 component observationally challeng- 
ing. Electron-impact excitation has been suggested as a 
means of directly probing the Ha in t he planet-forming 
regions of the disk (Bergin ct al. 2004). This mechanism 
requires stellar X-rays to create a distribution of photo- 
electrons that pump the molecules into excited electronic 
states. The H2 then relaxes, producing a characteris- 
tic far-UV cascade spectrum of discrete emission lines 
and quasi-conti nuous spectral fea t ures during the dis - 
sociation of IT (lAiello et al.l 119841: lAbgrah et aHll997f ). 
Ilngleby et al. (2009) proposed this mechanism to explain 
a portion of the low-resolution far-UV spectra from a 
sample of classical and weak-lined (low gas content) T 
Tauri stars observed with the Hubble Space Telescope- 
ACS and -STIS. Low-resolution and intermediate S/N 
in existing far-UV observations make it difficult to sep- 
arate the electron-impact H2 signal from Lya-pumped 
H2, atomic emission, CO, and the underlying accretion 
continuum. We lump these emissions together under the 
term "continuum" due to their mostly unresolved struc- 
ture in previous studies. Fortunately, the observational 
picture has improved dramatically with the installation 
of the Cosmic Origins Spectrograph (COS) on HST. 
HST-COS combines very low detector backgrounds with 
moderate spectral resolution (Aw w 17 km s _1 ). The low 
background permits measurements of the true continuum 
shape, including an assessment of the contribution from 
a hot accretion continuum. The resolution of COS also 
permits the identification and separation of the spectral 
components, enabling a more robust comparison with 
models of the processes that govern the planet-forming 
regions of CTTS disks. 

In this paper, we use COS to study in detail the far-UV 
continuum emission from these objects for the first time. 
We focus on two objects that we propose are prototypes 
for electron-impact H2 emission (V4046 Sgr) and contin- 
uum emission produced in the hot accretion shock (DF 
Tau). CO A - X band emission also contributes to the 
far-UV line and continuum spectrum in these systems. 
First results on CO emission will be presented in Paper 
II (K. France - in preparation). In §2, we describe the 
COS observations and data reduction. In §3, we describe 
empirical fits to the data and compare these fits with 
model H2 spectra. §4 describes our measurements of the 
collisionally excited H2 and accretion continua. We also 
discuss discrepancies between the new COS observations 
and predictions of our electron-impact H2 models. We 
then discuss in §5 the relevance of these processes to the 
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Fig. 1.— The far-UV spectra of DF Tau (black) and V4046 Sgr 
(green) in the 1140 — 1760 A bandpass, obtained by combining the 
G130M and G160M data and smoothed to 1 spectral resolution 
element for display purposes. Emission lines from hot gas (Si IV, 
C IV, He II, etc) and Lyo-pumpcd H2 dominate the spectrum in 
the far-UV bandpass. The flux is plotted in FEFU (= femto-erg 
flux unit, 1 FEFU = 1 X 10 -15 ergs cm -2 s" 1 A -1 ). The DF 
Tau spectra have been offset by +200 FEFU for clarity. The gray 
dashed line denotes the flux zero point in both observations. 



physical state of the disk and the local far-UV radiation 
field. Finally, we present a brief summary of our work in 
§6- 

2. TARGETS AND OBSERVATIONS 

DF Tau and V4046 Sgr are both binary pre-main se- 
quence systems with gas-rich disks. DF Tau is composed 
of two early M stars (0.68 and 0.51 M Q ; White & Ghez 
2001) with an age estimated to be between ~ 0.1 - 2 
Myr (see Ghez et al. 1997 for a discussion). DF Tau A 
is thought to be the stronger contributor to the far-UV 
output from the system with a mass accretion rate ap- 
pr oximately an order of magnitude greater than DF Tau 
B ([Herczeg et al.l 120061) . The DF Tau disk is observed 
approximately edge-on (i ~ 80 - 85°; Johns-Krull & 
Valenti 2001; Ardila et al. 2002). Emission from hot gas 
lines and Lya-pumped H2 have been studie d using exist- 
ing 5T-STIS echelle and FUSE spectra (jArdila et al.l 
120021 ; IHerczeg et al.l 120061) , and detailed analysis of the 
radiative transfer interplay between H I and H2 using 
COS will be presented by H. Yang (submitted). The 
0.1 - 4.5 keV luminosity of DF Tau is ~ 9 x 10 29 erg 
s _1 (jWalter fc Kuhil 119811). The outer disk has a dust 
mass ~ 130 M$ ([Andrews fc Williams 2005) and mm- 
wave CO has been detected ([Greaves! 120051). DF Tau 
show s evidence of a gas-rich inner disk ([Najita et al~l 
120031) . although a search for extend ed near-IR inner dust 
disk has returned a null detection ([Karr et~^[2OT0l) . 

V4046 Sgr is a close binary comprised of two mid-K 
stars (0.91 and 0.87 M ; Stempels & Gahm 2004) with 
a circumbinary disk at an intermediate inclination angle 
(i ~ 35°, Quast et al. 2000; Rodriguez et al. 2010). It 
has an age between 4-12 Myr, depending on its possi- 
ble membership in the P ic moving group (|Quast et al.1 
2000; IKastner et a l. 2008). It was s tudied as part of 
the IUE pre-main sequence star atlas ([Johns-Krull et all 
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TABLE 1 

HST-COS OBSERVING LOG. 



Object R. A. (J2000) Dec. (J2000) Date COS Mode T exp (s) 



DF Tau 04 h 27 m 02.81 s +25° 42' 22.1" 2010 January 11 G130M 4828 

DF Tau 04 h 27 m 02.81 s +25° 42' 22.1" 2010 January 11 G160M 5549 

V4046 Sgr 18 h 14 m 10.49 s -32° 47' 34.2" 2010 April 27 G130M 4504 

V4046 Sgr 18 h 14 m 10.49 s -32° 47' 34.2" 2010 April 27 G160M 5581 



2000). While FUSE data reveal H 2 emission from 
the Lyman and Werner levels, a search for mid-IR 
emi ssion from cooler Hg has returned only upper lim- 
its (jCarmona et al.ll2008ft . The 0.45 - 6.0 keV luminosity 
is estimated to be ^ 1 - 5 x 10 30 erg s _1 (Giinther et 
al. 2006; and §4.3). V4046 Sgr has a molecule-rich outer 
disk (M r s ~ HO Mq, M dust 40 ~ M e ; Rodriguez et 
al. 2010 L although the inner disk shows evidence of dust 
clearing ([Jensen fc Mathicu [1997D . The lack of near-IR 
excess categorizes V4046 Sgr as a "transitional" disk, al- 
though the ongoing accretion and far-UV molecular emis- 
sion suggest that a molecular gas reservoir is present in 
the inner disk. 

DF Tau and V4046 Sgr were observed with the medium 
resolution, far-UV (G130M and G160M) modes of COS 
on 11 January and 27 April 2010, respectively. They 
were each observed for a total of four orbits, with two or- 
bits per grating. In order to achieve continuous spectral 
coverage and minimize fixed pattern noise, observations 
in each grating were made at four central wavelength 
settings (A1291, 1300, 1309, and 1318 for G130M and 
A1589, 1600, 1611, and 1623 for G160M) in the default 
focal plane position (FP-POS = 3). This combination of 
grating settings covers the complete 1136 < A < 1796 A 
bandpass at a resolving power of R w 18,0003. Near- 
UV imaging target acquisitions were performed through 
the COS primary science aperture using MIRRORB. Ta- 
ble 1 provides a log of the COS observations acquired as 
part of this study. Table 2 summarizes relevant system 
parameters, including accretion diagnostics measured in 
this work. The data have been reprocessed with the 
COS calibration pipeline, CALCOSo v2.12, and com- 
bined wit h the custom IDL co addition procedure de- 
scribed bv iDanforth et all (|2010l ). Figure 1 displays the 
combined 1140 - 1760 A spectra of both systems. The 
flux units used on this and subsequent figures are FE- 
FUfS (1 FEFU = 1 x 10" 15 ergs cm" 2 s" 1 A" 1 ) . We 
use multi-epoch archival optical data spectra as support- 
ing observations of the spectral region around the B aimer 
break. 

3. ANALYSIS 

3.1. Emission from H2 and the Inner Accretion Disk 

The strongest features in the far-UV spectra of DF Tau 
and V4046 Sgr are the emission lines from hot gas and 



3 The COS LSF experiences a wavelength dependent non- 
Gaussian form due to mid-frequency wave-front errors produced 
by the polishing errors on the HST primary and secondary mir- 
rors; http: //www. stsci . edu/hst/cos/documents/lsrs/ 

4 We refer the reader to the cycle 18 

COS Instrument Handbook for more details: 

http : //www. stsci . edu/hst/cos/documents/handbooks/current/cos.coi^^^tmj, 

see §1.1.2 of the Cycle 17 COS Instrument Handbook 



photoexcited H2 (Figure 1). The hot gas lines (e.g., N V, 
Si IV, C IV, et c) are produced in th e magnetically active 
atmospheres (|Bouvier et al.l 120071) of the central stars 
and at (or near) the shock interface where the accretion 
strea m impacts the stellar surface (jGiinther fe Schmittl 
2008). The strong, discrete Lyman and Werner band 
H2 lines are mostly fluoresced by H2 pumping transi- 
tions that coincide with the broad Lya profiles in the 
these systems, altho ugh coincidences with other stellar 
lines also contribute (| Wilkinson et al. 2002; Fr ance et al.l 
l2007h . We see that the hot gas and photoexcited H2 lines 
are superimposed upon a mostly continuous, f aint un- 
derlying emission spectra. Ilnglebv et al.1 (|2009l ) argued 
that this emission is a combination of accretion lumi- 
nosity from near the hot accretion spots and emission 
from electron-impact excited H2. The electron-impact 
H2 can be further decomposed into discrete bound-bound 
emission and quasi-continuous emission from dissociative 
transitions to the ground electronic state, which is char- 
acteristic of the electron impact p rocess for electron dis- 
tributions with energies > 15 eV (lAiello et al.lll991lh In 
both systems, there is also a contribution from CO emis- 
sion which will be addressed in Paper II. In Figure 2, 
we show a breakdown of these components for V4046 
Sgr. The accretion continuum is represented as the blue 
dashed line, and a model H2 spectrum is shown super- 
imposed as a dotted red line. Because the vast number 
(several hundreds) of photoexcited H2 lines complicate 
the comparison of the broad band spectral characteris- 
tics of the data with the model, we isolate spectral re- 
gions that are free of known photoexcited H2 lines. These 
"continuum bins" have widths of ~ 0.6 - 1 A, and are 
shown as purple diamonds in Figure 2. 

The "continuum" spectra of the two systems are qual- 
itatively different. Figure 3 (left) shows the continuum 
spectrum of DF Tau, corrected for reddening (Ay = 0.6; 
H erczeg fc Hillenbrand 2008) using the interstellar curve 
of lCardelli et al.l (|1989t) . The far-UV data show a mostly 
linear rise with wavelength that we attribute to accretion 
emission. We show a linear fit to the data as the dashed 
red line and several optical spectra of DF Tau, represen- 
tative of the high (purple), intermediate (blue), and low 
(green) states of the object. We include a scaled ver- 
sion of the intermediate spectrum as the dashed black 
line to illustrate that the linear fit to the far-UV data 
could account for the observed continuous spectrum of 
DF Tau below the Balmer break (see §4.2 for a discus- 
sion of the spectral shape of this emission) . We attribute 
the dispersion of the far-UV points at A > 1500 A to 
contamination by unresolved H2 lines (both photo- and 
electron-excited), CO band emission, and the lower S/N 
data^af the red end of the far-UV bandpass where the 
'ective area decreases. The V4046 Sgr spectrum 
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TABLE 2 
System Parameters. 



Object 


Spectral Type 


Age 


Inclination 


logio {LxT 


Ref. 


logio (G cont /G ) h 


logio(£c/v) 


logio (M acc ) 






(Myr) 




(erg s" 1 ) 




(a = 1 AU) 


(erg s" 1 ) 


(M yr" 1 ) 


DF Tau 


MO. 5 + M3 


0.1 - 2 Myr 


> 80° 


29.9 


1,2,3,4 


6.8 


30.26 


-7.1 c 


V4046 Sgr 


K5 + K5 


4-12 Myr 


35° 


30.0 - 30.7 


5,6,7,8 


5.6 


29.25 


-7.9 C 



a For DF Tau, L x is the total 0.1 - 4.5 keV luminosity 
keV luminosity (see §4.3). 



v alter 



b Go is the strength of the average interstellar radiation field, evaluated over the 912 
Habing 1968). 

c Determined fr om the empir ical relationship betwee n C IV em i ssion and t he mass acc r etion rate, Eqn. (2 
Refere nces: (!) lArdila et al.l (20021) (2) I Johns-Krull fc Valentil OOOll), (31 IGhez et al.l JT997I '). (4) IWhite ' 
(2000), f6) IBtempels fc Ga hrn (2004), (7) IKastner et al.l ||200'8T) . (8) Rodrigue z et alj H20ldf ) 



For V4046 Sgr, L x is estimated for the total 0.45 - 6.0 
2000 A bandpass (Go = 1.6 X 10 



" erg cm 



oflJohns-Krull et al.l J2000D. 
Ghezi (120011 1, (5) IQuast et all 
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Fig. 2. — The binned continuum spectrum of V4046 Sgr (purple 
diamonds with error bars) is compared with a synthetic electron- 
impact H2 spectrum (blue squares) created for parameters from 
the literature [N(H 2 ) = 10 19 cm" 2 , T(H 2 ) = 4000K, E e = 20eV], 
The model has been offset by an underlying linear continuum, 
shown as the blue dashed line. The decomposition of the model 
into discrete and quasi-continuous emission is also shown (red dot- 
ted and pink solid lines, respectively). The discrete line emission 
has been shown at a decreased resolution for display purposes. We 
conclude that the observed continuum spectrum can be decom- 
posed into an accretion continuum (dashed blue line) and emission 
from H2. We note that the model overpredicts the observed emis- 
sion at A < 1300 A, as discussed in §3.3. 



includes emission with a shape qualitatively similar to 
that expected for electron-impact excited H2 (see Figure 
3, right). As we illustrated in Figure 2, this emission 
is superimposed upon an underlying linear continuum. 
We do not make an extin ction correction for V404 6 Sgr 
because Ay is neglig ible (IStempels fc Gahm|[200l . We 
consider DF Tau as an example of a CTTS with hot ac- 
cretion emission and V4046 Sgr as representative of a 
CTTS with electron-excited H2. In the following sub- 
sections, we describe the modeling procedure we used to 
develop synthetic H 2 spectra for a quantitative analysis 
of the electron-impact spectrum. The characteristics of 
each process are considered in more detail in §4. 

3.2. Synthetic Electron- Impact H 2 Spectra 

We created synthetic electron-impact H2 spectra 
for comparison with the CO S observations follow- 
ing the formalism described by Liu & Dalgarnol (|1996| ) 



and iWolven et al.l (|1997l) . Our spectral s ynthesis code 
includ es emission from photoexcited H2 (|France et al.1 
l2010bf ) , but we have decoupled the various emission com- 
ponents in order to study separately the discrete and 
quasi-continuous emission spectrum from collisionally ex- 
cited H2. We describe the calculation of electron-impact 
spectra here. The model includes emission from the B - , 
B - , B - , C - , D - , and D - X electronic tran- 
sitions, including Lyman band populations coming from 
E,F - B cascades. We explicitly include transitions to 
predissociating states and vibrational states that result 
in dissociation (v > 14, the vibrational continuum) and 
independently save the dissociation spectra. Transition 
probabilities and wavelengths fo r the B, B , C, and D 
were taken from the literature (|Abgrall et al.l fl993al fbl. 

I1994T ). the B and D parameters were calculated by 
combining the branching ratios of the B a nd D states 
with the appropriate energy level spacings (fKwok et al.1 
ll98aiGlass-MauieaiJl984tlHuber k Herzberelll979D . We 
computed transition probabilities for the E and F cas- 
cades using the Franck-Condon factors presented by iLinl 
(|1974h . The emission rate from level of nv' J' to level 
v" J" of the ground electronic state is given by 



S{n'v'J'^v"J") = R(n'v'J')[l-r]{n'v'J')] 



A 



(i) 

where R(n'v'J') is the excitation rate into the excited 
electronic level, i](n'v'J') is the efficiency for predis- 
sociation in the excited electronic state, and the ratio 
of t ransition probabilities (A nv j) is the branching ra- 
tio (|Liu fc Dalgarnolll996D . 

The excitation rate, R, is proportional to the column 
density in the ground electronic state, N v j, the cross- 
section connecting the ground state to the rovibrational 
level of the excit ed electronic sta te, a v j^ n ' v i j> , and the 
electron flux, F e (|Liu et al.lll99l) . 



R(nv'j') = F e N v ja v j^ n ' V 'j' 



(2) 



Given the uncertainties about the source of the free elec- 
trons and their geometry with respect to the spatial dis- 
tribution of the H 2 , we do not attempt to reproduce 
the absolute flux of the electron-impact emission in this 
work. Thus the electron flux can be assumed to be a 
scale factor. 

The cross-sections can be measured experimentally, 
but as we are not solving for the absolute flux of the 
electron-impact emission lines in the spectra of V4046 
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Fig. 3. — The far-UV and near-UV/optical accretion continuum in DF Tau (left) and V4046 Sgr (right). A linear fit to the far-UV 
spectra (red dashed line) is extrapolated to longer wavelengths, and compared with the "high", "intermediate", and "low" optical emission 
states (plotted in purple, blue, and green, respectively). The optical spectrum displayed as the black dashed line has been offset to show 
that a linear continuum can describe the accretion emission from the far-UV to the Balmer break. The extra emission above the dashed 
red line in the DF Tau data near 1600 A is most likely contamination from photo-excited CO and H2 emission, while for V4046 Sgr the 
extra emission above the linear fit is dominated by emission from electron-impact H2. 



Sgr, w e follow the procedure described by iLiu et all 
(1998) and calculate relative excitation rates described 
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by the equation 



J—^n'v' J' 



(3) 



where E e is the electron energy, uj v j is the degeneracy of 
the ground electronic state, and £l v j^. nv 'j> is the collision 
strength gi ven by the analy tic expression presented in 
Eqn (4a) of ILiu et all (post ). 

For simplicity, we assume a single population of col- 
lisionally excited gas. We input the ground state ther- 
mal temperature [T(H 2 )], the column density [N(H 2 )], 
and the electron energy. We created a grid of models 
in [T(H 2 ), A(iJ 2 ),_E e ]-space, including ten temperatures 
between 300 - 5000 K, five column densities in the range 
10 16 - 10 20 cm -2 , and ten electron distributions in the 
range 14 - 5000 eV. The lower limit at E e > 14eV is set 
by the minimum energy needed to populate the C, E, 
and F levels in our model. Numeric errors occur when 
these states are not appreciably populated. Physically, 
the absolute cross-sections for electron excitation into the 
Lyman and Werner bands are very low at E e < 14eV, 
with oh2 % 3 x 10~ 18 cm 2 at E e < 12eV (compared to 



&H2 > 

1998). 



2 x 10 



-16 



cm 2 at 20 < E e < 200eV; Liu et al. 



3.3. Comparing the H 2 Models with the V4046 Sgr 
Spectra 

The model spectra computed on a 0.01 A grid are 
convolved with the appropriate COS linespread func- 
tion (LSF) for comparison with our spectroscopic ob- 
servations. The electron-impact spectrum of H2 con- 
sists of both discrete and quasi-continuous emission 
features. The continuum features are strongest at 
A > 1350 A (jAbgrall et all 119971) . while the electron- 
impact spectral energy distribution is mostly discrete at 
the blue end of the COS bandpass (A < 1300 A). Figure 
2 shows the qualitative agreement of the V4046 Sgr data 
with the H2 model, and we anticipate that a broad-band 
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Fig. 4. — A comparison of the full-resolution short- wavelength 
observations of V4046 Sgr (black) with the predicted discrete line 
emission from an unattenuated electron-impact H2 model (blue 
dashed line; see §3.3). Despite the agreement at A > 1400 A, many 
of the predicted features are not observed. Atomic species (e.g., 
N I, Si III, and N V) are labeled in purple, and emissions from 
photoexcited H2 are labeled in orange. The photoexcited lines are 
listed in Table 2. 



X 2 minimization approach can be used to estimate the 
H2 parameters of the system. 

In order to explore what additional processes should 
be included in the models prior to fitting the data, we 
compare the broad spectral colors with the intrinsic color 
ratio (CR) used to characterize electron penetration into 
gas giant atmosp heres, CR = 7(1550 - 1620 A)/I(1230 - 
1300 A) - 1.1 (|Yung et al.1 H98l IGustin et al.ll2006l) . 
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TABLE 3 

V4046 Sgr photoexcited H2 emission lines in Figures 4 and 6. 




Line ID a 


^oba 

(A) 


Arest 

(A) 


Pumping Transition 13 


^pump 

(A) 


C (1-4) P(10) 
B (4-1) P(14) 
C (1-4) P(ll) 
B (1-2) R(3) 
B (3 - 0) P(17) 
C (1-5) R(3) 
C (1 - 5) Q(3) 


1193.37 
1197.15 
1198.00 
1202.44 
1203.24 
1205.53 
1208.90 


1193.34 
1197.14 
1197.99 
1202.45 
1203.24 
1205.53 
1208.93 


(1 - 5) R(8) 
(4 - 2) R(12) 

(1 - 5) R(9) 

(1 - 2) P(5) 
(3 - 1) R(15) 

(1 - 5) P(5) 
(1 - 1) Q(3) c 


1214.62 
1213.68 
1217.00 
1216.07 
1214.47 
1216.99 
1031.86 


C (3-6) P(17) 
C (4 - 8) R(12) 
B (4 - 4) R(3) 
C (1 - 6) Q(3) 
C (4-8) P(8) 
B (1 - 3) R(3) 


1249.99 
1252.05 
1253.63 
1254.11 
1255.55 
1257.81 


1250.00 
1252.10 
1253.67 
1254.11 
1255.56 
1257.83 


(3 - 5) P(17) 
(4 - 7) R(12) d 
(4 - 3) P(5) 

(1 - 1) Q(3) c 
(4 - 7) P(8) 
(1 - 2) P(5) 


1216.74 
1220.47 
1214.78 
1031.86 
1219.95 
1216.07 



a B denotes a Lyman band (B x Sj — X transition, C denotes a Werner band (C 1 T1 U — -^^g") transition. 
b H2 lines excited by the broad stellar H I Lycr emission line, unless otherwise noted. 
c Pumped through a coincidence with the stellar O VI A1032 A resonance line. 

d The observed line may contain additional H2 fluorescence excited through coincidence with H I Ly/3 and/or Si III A1206.50 A. 



a 1.0 



o 
Z 



0.5 



0.0 




H, Transmission 
CH 4 Transmission 
e-lmpact H, 
Transmitted 
H 2 + CH4 
log l0 N(H,) = 19.00 
T(H,) = 4000 K 



1.0 



0.5 



1100 1200 



1300 1400 1500 
Wavelength (A) 



1600 1700 



Fig. 5. — (top) Transmission functions of methane and optically 
thick H2 through the disk of V4046 Sgr. (bottom) The synthetic 
electron-impact H2 spectrum is plotted in blue, the synthetic spec- 
trum attenuated by H2 self-absorption is shown in green, and the 
model attenuated by H2 + CH4 is shown in magenta. The blue, 
green, and magenta curves are the same at A > 1400 A, showing 
that molecular absorption most strongly impacts the observations 
in the G130M bandpass. 



where 



J(AA) = 



/■A; 



F H2 d\ 



(4) 



and is the electron-impact H 2 flux. Confusion from 
the overlap with hundreds of strong photoexcited H2 
emission lines in the CTTS spectra forced us to eval- 
uate the data in m 1 A bins, as described in §3.1. Using 
this binning, we find an intrinsic model color ratio of 
CRmod ~ 1-33. The observed color ratio in V4046 Sgr 

6 However, when our models are evaluated over the entire "blue" 
and "red" bandpasses for reasonable input parameters, we find 
that the total integrated flux is actually slightly greater is the sub- 
1400 A band [7(1140 - 1400 A) > 7(1400 - 1680 A)], which is a 
consequence of the large contribution from discrete features. This 



is CR bs = 2.1 ± 0.2, indicating that the observations 
are significantly redder than their intrinsic colors. Jo- 
vian color ratios in this range are typically int erpreted as 
internal reddening in th e upper atmosphere (iDols et aLl 
211(111: IGustinlt^[200l . 

The model predicts the spectrum of the discrete lines 
that have been observed in short -wavelength far-UV 
spect ra of the gas giant planets (|Gustin et al.l 12004 
2009). As a first-order check on the discrete lines, 
we compare in Figure 4 a section of the full-resolution 
COS G130M data with an electron-impact H 2 model. 
The model, overplotted as the dashed blue line, as- 
sumes parameter values typical of previous electron- 
impact studies of CTTS disks [iV(H 2 ) = 10 19 cm- 2 , 
T(H 2 ) = 4000K, E e = 20eV], although this electron 
distributi on has a sligh t ly hig her energy than those fa- 
vored by llngleby et al.1 (|2009f ). We discuss the electron 
energy distribution in §4.3. Figure 4 shows that many 
discrete features predicted by the model are not seen in 
the data. Because of this disagreement and the redden- 
ing of the broad-band colors, we conclude that calcu- 
lations of electron-impact spectra taken from the liter- 
ature are alone inadequate for direct comparison with 
the COS observations. In the following subsections, we 
therefore develop a more detailed physical picture that 
produces agreement with both constraints provided by 
the observed H2 spectrum. Table 3 presents a list of the 
photoexcited H 2 lines identified in Figure 4. 

3.3.1. H2 Self-Absorption 

A strong motivation for studying the electron-impact 
excited H2 spectrum in protoplanetary disks is that the 
more easily observable Lyq-pumpe d H2 emission is likel y 
produced on the surface of the disk ([Herczeg et al.l f2006) , 
whereas the X-rays needed to produce the non-thermal 
electron populat ion can penetrate deeper into the plane t 
forming region (jBergin et al.l 12004 unglebv et al.ll2009|L 
Thus, the electron-impact excited gas may be a use- 
qualitative result is true for our model spectra over a wide range 
of physical parameters, that is, the "blue" and "red" fluxes are 
approximately equal. 
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ful probe of the molecules that are in the actively 
planet-forming regions of the disk. The role of radia- 
tive transfer in determining the resulting UV spectrum 
measured at Earth has been consid ered for radiative ex- 
citation of H2 in the disk surface (Hcrczcg 
but not for the collisionally excited H2 in the inner 
disk. We first consider reddening by dust. We dis- 
count strong internal reddening because grains in the 
disk have likely grown ap preciably during the few Myr 
age ty pical of CTTS disks (jFurlan et aLll2005t lApai et al.l 
2005), and there is evidence that large grains produce 
a rela tively "gray" extinction ac r oss the far-UV band- 
pass (jFitzpatrick fc Massal 119881: iCardelli et all [1989). 
Grain growth and settling has also been shown to sig- 
nificantly reduce the ove rall opacity to UV p hotons in 
disks older than ~ 1 Myr (|Vasvunin et al.ll20l6T ). Also, if 
these inner disk grains play a large role, they must also 
extinguish stellar emission in edge-on systems, which is 
clearly not the case (e.g., the strong continuum and hot 
gas lines are observed in DF Tau). 

An effect that must be considered is self- absorption 
by H2 in the outer layers of the disk. At first con- 
sideration, it seems plausible that H2 photons emitted 
in the disk plane are re-absorbed before they can es- 
cape. This would naturally explain both the red CR and 
the absence of discrete emission features because self- 
absorption would remove individual lines from the re- 
sultant emission spectrum, thereby reddening the model 
spectrum. As cool H2 [T(ti2) < 700 K] has apprecia- 
ble optical depth only in the first few rotational levels 
(J ~ 6), this gas will not produce significant absorption 
out of excited vibrational levels. This cool H2 resem- 
bles the well-studied spectrum of H2 in interstellar clouds 
(e.g., Rachford et al. 2002; Burgh et al. 2007) and can- 
not provide appreciable opacity in the COS bandpass 
(A > 1150 A). A population of hot H 2 [T(H 2 ) > 2000 
K] would be required. We have tested this hypothe- 
sis by creating a tunable absorbing "screen" of H 2 that 
will attenuate the radiation produced by the electron- 
impact model. We created H2 absorption spectra for 
2000 < T(H 2 ) < 5000 K by calculating the Boltzmann 
distribution for the ground electr onic state and used the 
H 2 optical depth templates from iMcCan dliss (2003|) to 
compute transmission functions as a function of temper- 
ature and column density. We included all rovibrational 
states in the range (0 < J < 14) and (0 < v < 3). Even 
larger temperatures would be required to thermally pop- 
ulate higher states, but we do not consider such temper- 
atures because they are signifi cantly above the » 4 500 K 
H 2 dissociation temperature (jLepp &: Shulll IT983T ) . Fig- 
ure 5 shows the transmission function of H2 in red at the 
top, while the intrinsic (blue) and transmitted (green) 
H2 emission spectra are shown at the bottom. 

In general, we find that H2 self- absorption is in- 
sufficient to account for the color and discrete fea- 
ture discrepancies described in §3.3. Our results for 
H2 self- absorption of discrete lines are shown as illus- 
trative examples in Figure 6 for typical temperatures 
[T(H 2 ) = 4000K] and column d e nsitie s [N(R 2 ) = 10 19 
cm~ 2 ] obtained by llngleby et al.l ()2009l ). The top panels 
of Figure 6 show the computed H 2 transmission spectrum 
in red, the unattenuated electron-impact model in blue, 
and the transmitted H 2 spectrum in green. We find that 



while several of the transitions are indeed completely ab- 
sorbed, many of the strongest transitions are transmit- 
ted. It is these emission lines that we would expect to 
observe. The lower panels of Figure 6 show that these 
emissions are not observed, even though the COS sensi- 
tivity curve peaks in the A w 1250 A window displayed in 
Figures 4 and 6. The primary reason that self- absorption 
alone cannot explain the discrepancy is that the H2 pop- 
ulation created by electron-impact is highly non-thermal, 
whereas upper layers of the disk are mostly described by 
a warm thermal distribution (2000 < T(H 2 ) < 3000 K; 
Herczeg et al. 2006) with deviations caused by UV pho- 
toexcitation. We conclude that while H 2 self-absorption 
must be included, it is not the dominant factor in de- 
termining the shape of the emitted electron-impact H2 
spectrum in V4046 Sgr. 

3.3.2. Hydrocarbon Absorption 

In the previous section, we discussed the contribution 
of dust and H 2 self-absorption to the spectral characteris- 
tics of the observed spectrum of V4046 Sgr; in particular, 
that the long-wavelength (A > 1450 A) continuum could 
be described by emission from H2 excited by non-thermal 
electrons, but there is a lack of discrete emission lines 
at shorter wavelengths. We have determined that nei- 
ther dust nor H 2 self- absorption can attenuate synthetic 
electron-impact H2 spectra to the extent observed. 

We now consider absorption by more complex 
molecules, CO, H2O, and hydrocarbons. We first as- 
sume that the electron-excited H 2 is located in the disk 
and that the H2 emission must pass through a "screen" 
of these more complex molecules to reach an external ob- 
server. We require that the attenuating screen removes 
the short-wavelength H2 emission, but not the flux at 
longer wavelengths. While CO is observed in UV spec- 
tra of CTTS disks (e.g., Paper II), the A - X spectrum 
that is the most prominent CO feature at far-UV wave- 
lengths has a characteristic band structure. Even for 
very high temperatures, the interband spectral regions 
have low opacity. As a simple quantitative test for the 
importance of CO absorption within the disk, we cre- 
ated a template CO optical depth spectrum for column 
density iV(CO) = 10 17 cm" 2 at T = 1000 K, with spec- 
tral binning of 0.05 A. Even for this warm, high col- 
umn density medium, we find that only approximately 
7% of the spectral bins have tco > 0.1. Thus, there 
are ample low-opacity windows through which the un- 
derlying H 2 spectrum can escape. A similar argument 
can be made for the spectrum of water. Laboratory 
studies have identified individual vibrational modes of 
the H2O Rydberg ser ies in the 6-11 eV energy re- 
gion (Mo ta et al.| [2005). The far-UV cross-sectional spec- 
trum is characterized by strong bands in the 9.9-10.8 eV 
(1250 > A > 1150 A) region where we require substantial 
attenuation of the H2 electron-impact spectrum. How- 
ever, because the interband spectral regions have cross- 
sections of <7ff20 2 x 10~ 18 cm~ 2 , very large column 
densities of water would be required for attenuation of 
the electron-impact H 2 emi ssion which is thou ght to be 
produced in the inner disk ([Bergin et al.ll200~4l ). The re- 
quired column densities needed to reach TH20 ~ 1 are 
larger than those found in observational and model stud- 
ies of H 2 in the mid-IR spectra of protoplanetary disk 



France, Yang, & Linsky 




1240 1245 1250 1255 1260 

Wavelength (A) 



Fig. 6. — A combination of Figures 4 and 5 illustrating the effects of H2 and CH4 absorption at full spectral resolution. The V4046 Sgr 
spectra are plotted in black. The red H2 transmission curves show that even at high temperatures and relatively high column densities 
[T(H2) = 4000 K, N(H2) = 10 19 cm~ 2 ] there are significant windows through which H2 emission can escape in the COS G130M bandpass 
(shown in green). By contrast, more continuous absorption by hydrocarbons can effectively remove the short wavelength flux from H2 
(shown in magenta) while allowing the electron-impact excited H2 emission observed in the COS G160M bandpass to escape (see Figure 
5). Assuming that all of the hydrocarbon absorption can be attributed to CH4, we find that JV(CH4) > 6 X 10 16 cm -2 is required to 
suppress the discrete line emission to the observed level. 
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systems (jCarr fc Naiital [20081 iBethell fc Bergml l2009h . 
Finally, H2O h as a prominent absorption band in th e 
7 - 8 eV range (|Diercksen et al.lll982t iMota et al.ll2005D . 
which is coincident with the broad 1575, 1608 A spec- 
tral features characteristic of electron-impact excitation 
of H2. This dissociative H2O absorption feature has ab- 
solute cross-sections that are comparable to or stronger 
than those in the short-wavelength windows described 
above. Thus water column densities high enough to at- 
tenuate the bluest wavelengths in our data would also 
eliminate the flux at 1600 A, which is not observed. 

We suggest that a more likely scenario involves the ab- 
sorption of H2 emission by hydrocarbons in the planet- 
forming disk. This scenario is analogous to what has 
long been observed in far-UV auroral spectra of gas 
giant planets: a mostly electron-impact H2 spectrum 
that is "reddened" by hydrocarbon absorption, mostly 
from methane, ethane, and acetylene (CH4, C2H6, and 
C 2 H 2 ; Clarke et al. 1980; Yung et al. 1982). In par- 
ticular, CH4 and C2H6 have absorption cross-sections 
that meet the required criteria: (i) strong (er > 10~ 17 
cm -2 ), quasi-continuous absorption below A = 1400 A, 
and (ii) negligible effect on the output spectrum (er < 
10~ 18 cm~ 2 ) at the long-wavelength end of the COS 
bandpass. In the Jovian atmosphere, methane is often 
assumed to be the most important species reg ulating the 
broad -band far-UV spectral sig nature of H2 (|Kim et al.1 
Il997t IGrodent fc GerardI [200l . We compute the total 
methane opacity required to match the observed flux lim- 
its from the non-detection of discrete H2 electron-impact 
features in the COS data at A < 1300 A. We first scaled 
the 1575 - 1608 A model emission spectrum to match 
the observed flux. We then included self-absorption fol- 
lowing the procedure described above, and then applied 
a methane transmission function to the data. We find 
that iV(CH4) >6x 10 16 cm -2 is sufficient to attenuate 
the short-wavelength discrete features within the errors 
of the observed V4046 Sgr data. The relative abundances 
of the hydrocarbons in the disk are unknown, and a sim- 
ilar exercise assuming a purely ethane screen requires 
A(C 2 H 6 ) > 4 x 10 16 cm A The difference is due to the 
somewhat larger cross-sections of ethane in the 1140 - 
1300 A bandpass. 

We caution the reader that our assumption of pure 
methane absorption in the disk was chosen by analogy to 
planetary atmospheres. Although protoplanetary disks 
are the building blocks of these atmospheres, the anal- 
ogy is uncertain. We do not have a strong observa- 
tional constraint on the choice of methane over ethane, 
although the former is almost always found to be the 
dominant source of far-UV attenuation in the Jovian at- 
mosphere (Dols et al. 2000). A combination of hydrocar- 
bons is undoubtedly a more likely model, and there may 
also be additional species that we have not considered 
in the present analysis. Whatever the exact makeup of 
the disk, some sort of absorbing screen is required if the 
usual assumption that H2 electron-impact emission is a 
combination of bound-bound and dissociative transitions 
is valid in these disks. Modeling the far-UV radiative 
transfer in these disks is required to obtain a more quan- 
titative understanding of the spectra we now routinely 
observe with HST-COS. 
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Fig. 7. — The best-fit model of the electron-impact excited H2 
emission in the spectrum of V4046 Sgr. Due to the attenuation 
of the short-wavelength discrete emission, we restrict the fits to 
A > 1400 A. The best fit model parameters are T{H 2 ) = 300Q+goo° 



N(H 2 ) 



10 8 cm , and electron energies in the range 



K, 

E e ~ 50 - 100 eV. The H2 column density is only weakly con- 
strained by this analysis. We also show the low energy model 
(E e = 14 eV) for the best fit temperature and column density. 
The underprediction of the model flux near 1600 A is due to not 
including CO emission from the disk. 



4. DISCUSSION 

4.1. Electron- Impact H2 Emission in V4O46 Sgr 

In §3.3, we described how the A < 1400 A discrete 
emission lines from the electron-impact H2 excitation 
cascade are absorbed and/or scattered out of the line- 
of-sight. In order to measure the properties of the colli- 
sionally excited H2 spectrum, we performed a x 2 mini- 
mization analysis to compare the binned V4046 Sgr data 
at A > 1400 A with the grid of electron-impact H2 mod- 
els described in §3.2. The free parameters are [JV(H2), 
T(H2), E e ] and a linear fit to the underlying continuum, 
which we will discuss in the following subsection. We find 
that the electron-impact excited H2 has a temperature of 
T{H 2 ) = 30001^8° K, a column density of 7Y(H 2 ) ~ 10 18 
cm~ 2 , an electron energy of E e ~ 50 - 100 eV, but 
the only firm constraint is on the temperature. By re- 
stricting the fits to A > 1400 A, we avoid the 1200 - 
1300 A region where our models predict a dependence 
on A(H 2 ). We find a minimum near E e ~ 50 - 100 eV, 
but higher energies are generally consistent with the data 
for 7V(H2) = 10 16 ~ 19 cm~ 2 . In Figure 7, we show repre- 
sentative fits to the A > 1400 A data. We find that the 
best fit parameters not only display the lowest reduced 
X 2 , but they also show reasonable fits to the eye. For 
comparison, we also display a low energy case (E e = 14 
eV). 

A major complication to all of these fits is the spectral 
overlap between the electron-impact emission spectrum 
of H 2 and CO A - X band emission. With our broad 
binned approach, we can avoid strong discrete lines from 
the stellar atmosphere and accretion shock, as well as 
photoexcited H2 . The CO emission is distributed in wide 
vibrational band structures with non-zero intraband flux. 
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Using the binned continuum bands described above, we 
necessarily include CO emission. We will present models 
of this emission in Paper II. For now, we note that the 
strongest contaminant to the pure H2 electron-impact 
spectrum is the CO (0 - 1) A1597 A band that fills in 
the spectral valley between the H2 emission peaks at 1575 
and 1608 A. This excess emission tends to push the fits 
to higher temperatures, where the dissociation contin- 
uum broadens, or to very low electron energies (E e < 14 
eV), where dissociation is less important. This may ex- 
plain some of the discrepancies between our results and 
the best-fit parameters \T(H 2 ) ~ 5000 K; E e ~ 12 eV] 
presented bv llnglebv et all (|2009t) . 

4.2. Accretion Continuum 

Figure 3 shows the spectra of our targets from the 
far-UV to the optical. The black dotted optical spec- 
trum has been scaled from an archival observation to 
show that the linear fit to the far-UV data alone (dashed 
red line) can plausibly also describe the shape of the 
A < 3600 A Balmer continuum used to measure the ac- 
cretion i^rtes_in_low-m^ 

disks (ICalvet fe Gullbrinel H99l IHerczeg et all f2009h . 
We fit the far-UV spectra with both a power-law and 
a straight line. For both spectra a simple linear fit of 
the form F\ = m\ + b produced the best match to the 
data. We find [m,b] = [0.0393, -44.0250] and [0.0026, 
0.0918] for DF Tau and V4046 Sgr, respectively, where 
the units of m and b are FEFU A -1 and FEFU, where 1 
FEFU = 1 x 10~ 15 erg cm~ 2 s" 1 A -1 . The linear slope 
for DF Tau is ~ 15 times that found for V4046 Sgr. 
While DF Tau is known to have a much larger accretion 
rate (2 < M acc < 1300 x 10~ 9 M Q yr" 1 ; Hartigan et al. 
1995; Herczeg et al. 2006; Herczeg & Hillenbrand 2008) 
than V4046 Sgr (M acc - 3 x 10~ n M Q yr^ 1 ; Giinther 
and Schmitt 2007) a direct comparison between accre- 
tion rates from the literature and the far-UV accretion 
continuum slopes is unreliable given the different mea- 
surement methods and temporal variability. 

We can estimate the mass accretion rates directly from 
the COS observations using empirical scaling relations 
between the C IV luminosity (Lciv) an d M acc . The 
largest systematic source of uncertainty to this approach 
is the relative contribution from the stellar transition 
region emission t o the observed C I V pro file. We use 
equation 2 from Uohns-Krull et ahl (|2000t ) to estimate 
the mass accretion rate in these systems at the time 
of our continuum observations, assuming that the C IV 
surface flux in excess of a saturated stellar component 

(Fciv 
shock. 

We fit the C IV profiles with a multi-component line- 
profile employing the appropriate COS LSF. Assuming 
distances of 140 and 70 pc for DF Tau and V4046 Sgr, 
respectively, we find total L C iv = 1.81 ± 0.09 x 10 30 erg 
s _1 and 1.76 ± 0.04 x 10 29 erg s _1 . The error bars are 
the measurement errors, which are very small for these 
bright lines. These are the time-averaged luminosities 
over the total G160M observing time. The C IV fluxes 
were found to be slowly varying over the » 3 hours of 
G160M exposures. We measure C IV count rate changes 
of +7.0 and -13.6% for DF Tau and V4046 Sgr from the 
COS microchannel plate two-dimensional spectrograms, 



> 10 6 erg cm 2 s x ) is produced at the accretion 



evaluated in 50s time intervals. 

These C IV luminosities can be explained by mass ac- 
cretion rates of 7.8 x 10^ 8 M Q yr _1 and 1.3 x 10~ 8 M© 
y r -i £ or pjp Tau and V4046 Sgr, respectively. While our 
accretion rate fo r DF Tau is only slightly large r than the 
range quoted bv IHerczeg fe Hillenbrandl ((2008), our mea- 
surement for V4 046 Sgr is considera b ly lar ger than the 
value quoted by IGiinther fe Schmittl (120071) . Th e large 
discrepancy with the IGiinther fe Schmittl (|2007f ) accre- 
tion rate may be attributable to a larger contribution 
from transition region C IV in V4046 Sgr, however it 
seems unlikely that values of the mass accretion rate as 
low as Macc ~ 3 x 10~ n Af© yr -1 are consistent with 
our observations of V4046. This is not entirely surpris- 
ing as CTTS mass accretion rates based on X-ray obser- 
vations typically differ from optical/UV determinations 
by more than an order of magnitude due to absorption 
in the high-density regions where shock X-rays are pro- 
duced (see e.g., Sacco et al. 2010 and references therein). 
The relevant accretion measurements are given in Table 
2. 

We suggest that future observations for which far-UV, 
near-UV, and optical observations are obtained in close 
time proximity are required to measure accurately both 
the far-UV continuum and the accretion rate from the 
Balmer continuum. These observations will tie together 
the observed far-UV accretion continuum in these objects 
with the well-calibrated accretion rate diagnostic. 

4.3. Outstanding Issues Regarding Electron- Impact 
Excitation of Hi 

There is some uncertainty in the assumption that 
electron-excited H2 is the dominant molecular contin- 
uum emission mechanism for V4046 Sgr and other young 
circumstellar disks. Our models and analysis support a 
picture in which intermediate energy electrons excite H2 
in the inner disk, which then radiates and is observed 
superimposed upon an underlying continuum from the 
accretion shock. Since X-ray s are invoked to create the 
photoelectrons (fDrainc 1978J) needed to collisionally ex- 
cite the H2, we must ask if X-rays are sufficient for 
this mechanism to operate. We can measure the to- 
tal flux from electron-excited H2 in V4046 Sgr, F112 ~ 
8 ± 4 x 10 -13 erg cm -2 s^ 1 (approximately 2 x the 
integrated 1400 - 1660 A H2 flux shown in Figure 7 in 
order to account for the implied A < 1400 A emission, 
§3.3). If we assume an inner disk electron-impact H2 
emitting region (0.2 < ani < 0.5 AU; see §5.2 for addi- 
tional discussion) , the collisionally excited H2 luminosity 
is L c °f ~ 9 ± 5 x 10 28 erg s" 1 , or ~ 6 ± 3 x 10 39 
photons s _1 for a fiducial H2 wavelength of 1400 A. 

High resolution Chandra spectra of V4046 Sgr show 
that the accretion shock flux is appr oximately half 
of the observed value in T W Hya (Giinther et al.l 
l200i IGiidel fe Telieschl [2007h . Scaling from the ob- 
served _0 ; 45_j^_6 ; 0keVX-ray luminosity (L x ) of TW 
Hya (Ka stner et al.l I2002f) . we estimate the total X- 
ray luminosity for V4046 Sgr is Lx ~ 1 - 5 x 10 30 
erg s _1 , or 0.6 - 3 x 10 39 photons s^ 1 for a fidu- 
cial X-ray energy of 1 keV. 1 keV is the approxi- 
mate peak energy o f the V4046 Sgr X-ray spectral 
energy distribution (| Giinther fc Schmittl l2007f ). The 
grain photoelectric yields for the sub-micron sized grains 
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that are present at the inner edge of transitional dust 
disks ([Eisner et al.ll2006f) are of order un ity for photons 
with E>1 keV ([Weingartner et al.ll2006[ ). Thus, we find 
that V4046 Sgr produces sufficient X-ray luminosity to 
create enough photoelectrons to excite the observed H2 
emission from the inner disk. 

While this electron-impact H2 interpretation is en- 
ergetically plausible for the observed emission, it may 
not be adequate to describe the full far-UV emission 
spectrum of V4046 Sgr. Absorption by a hydrocarbon 
screen appears to reproduce both the observed absorp- 
tion of the short-wavelength discrete lines and the un- 
obscured H 2 emission at A > 1400 A, but verification of 
this result requires additional study. A more specula- 
tive scenario involves a very high column density, high 
temperature self-absorbing screen of H2 that could pro- 
duce nearly continuous absorption in the COS G130M 
band while not absorbing the long wavelength emission 
arising from dissociative transitions. Another possible 
scenario is a "non-traditional" electron-impact H2 emis- 
sion process in which dissociative yields of order unity 
are achieved. This process would produce the quasi- 
continuous long wavelength emission without most of 
the energy being distributed among the discrete states 
of the Lyman and Werner bands. As the canonical H2 
dissociation fraction from the L yman and Werner bands 
is ~ 0-15 dStecher fc Williams! Tl967t IShull fc Beck^ithl 
Il982t lAbgrall et al.lll997l) . large dissociative yields would 
require substantially fine-tuned molecular and electron 
distributions. It seems unlikely that either of these sce- 
narios are physically possible. The other prominent long- 
wavelength emission source in these disks is CO, although 
CO alone probably cannot account for the observed far- 
UV molecular spectra of CTTS disks. 

We point out that a spectrograph such as COS is es- 
sential for a comprehensive analysis of such systems. The 
low background and high-sensitivity are both important 
because the continua in these systems have characteristic 
flux levels of a few x 10~ 15 erg cm~ 2 s _1 A -1 or less. 
Moderate spectral resolution is also required for the de- 
tection of individual discrete lines predicted by electron- 
impact models. Moreover, when protoplanetary systems 
are observed without heavy internal reddening, spectral 
coverage at A < 1300 A is important because this re- 
gion probes the discrete line emission from the Werner 
bands of H2 whose spectrum is a more sensitive diag- 
nostic of the electron energy distribution. Without this 
wavelength region, spectral fits to the data will not ac- 
count for the Werner lines, biasing the results toward 
very low energy electron distributions. As the threshold 
energies f or the Lyman and Werner bands are 11.37 and 
12.41 eV (jAiello et al.lll984T ) . respectively, Lyman band 
continuum emission in the absence of discrete Werner 
band lines requires a highly-tuned electron distribution. 
Because Lya productio n from the disso ciation of H2 
drops sharply at 14.7eV ([Aiello et al.ll991| ). discrete lines 
should be even more important relative to the continuum 
a t very low elect r on ene rgies . 

iHerczeg et al.l (|2006i) presented high-quality HST- 
STIS E140M observations of DF Tau. While the spec- 
tral resolution of these data is sufficient to analyze the 
brighter Lya pumped H2 emission, the accretion and H2 
continua are too confused by the instrumental noise floor 




1100 1200 1300 1400 1500 1600 1700 
Wavelength (A) 

Fig. 8. — Comparison of the spectral shape of the far-UV accre- 
tion shock emission in DF Tau with photospheric emission from the 
Herbig Ae star AB Aur. The linear fit to the DF Tau continuum is 
shown as the red dashed line. The strong continuum emission from 
DF Tau is an effective radiation field that is ~ 3 times hotter than 
the radiation from the primary star, creating a source of far-UV 
continuum photons in the planet forming region of the disk. The 
enhanced flux (above the linear fit) in DF Tau near 1600 A is most 
likely CO emission. 



to ma ke a clear separation. Conversely, llnglebv et al.l 
(2009) presented a large sample of CTTS targets using 
the ACS-SBC, but the low spectral resolution makes the 
detection of discrete molecular features impossible. Nei- 
ther of these instruments has sufficient throughput in the 
A < 1300 A band to measure the Werner emission lines 
of H 2 . 

5. PROBING THE FAR-UV RADIATION ENVIRONMENT 
AND MOLECULAR GAS OF INNER CTTS DISKS 

5.1. The 912 - 2000 A Radiation Field and Molecular 
Dissociation 

The high sensitivity and spectral resolution of COS al- 
lows us to isolate the accretion shock continuum in pre- 
main sequence circumstellar disk systems. This higher 
temperature continuum radiation fi eld is similar to that 
discus sed for longer wavelengths bv lCalvet fc Gullbringj 
(1998) and the multi-comp o nent fit to the UV spectrum 
of TW Hya bv lCosta etaLl (|2000f ). In Figure 8, we show 
the continuum spectrum of DF Tau, whose central ob- 
ject is an M2 binary system (see §2 ), with a photospheric 
temperature of TeU ~ 3500 K ([Hartigan et al.l [19951 : 
Uohns-Krull et al. 2000). The effective temperature of 
the accretion emission in DF Tau is « 10 4 K, which 
creates a local radiation field in the planet-forming re- 
gion of the disk with a shape similar to that of a Herbig 
Ae star. To illustrate this point, we compare the scaled 
iJST-STIS spectrum of AB Aur in Figure 8 (green) with 
the binned DF Tau continuum spectrum. AB Aur and 
DF Tau have similar ages. AB Aur is known to be sur- 
rounded by a molecule-rich disk (H2 and CO; Roberge 
et al. 2001), and t here is evidence that plane ts are form- 
ing in the system (fOppcnhcimc r~et al.ll2008i ). A critical 
difference between DF Tau and AB Aur (spectral type 
AOVe, T eff w 10 4 K) is the temperature of the UV radi- 
ation field. 
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The accretion radiation detected down to A = 1147 A 
with high confidence provides strong evidence for contin- 
uum emission at wavelengths down to the Lyman limifl 
The shape and intensity of this far-UV radiation field 
have a significa nt impact on the evolutio n of a proto- 
planetary disk ([van Dishoeck et al.l 120081 ) . In particu- 
lar, dissociation of the most abundant species is con- 
trolled by photons in this critical 912 - 1150 A wave- 
length r egion. H 2 is photodissocia ted through a two-step 
process (jStecher fc Williamslll967l ) starting with absorp- 
tion into the Lyman bands at A < 1108 A ()Carruthers! 
Il970t) . extending to longer wavelengths for rovibra- 
tionally excited molecules. The first strongly predis- 
sociating CO absorption system is the E - X band 
at A < 1076 A (|van Dishoeck &; Blackll988D . with shorter 
wavelength far-UV transitions dominating the dissoci- 
ation as radiative transfer effects become important 
deeper into the disk. The dissociation of H 2 and CO 
molecules produce species that can catalyze the chem- 
istry throughout the inner disk, including neutral H, vi- 
brationally excited H2 (H|), neutral and ionized C, and 
O. This suggests that the boundary regions where the 
far-UV continuum penetrates the higher density molec- 
ular disk could be considered photodissociation regions 
(PDRs). Far-UV continuum emission will contribute to 
the dissociation of water (as discussed in §3.3.2), which 
is relevant for both giant planet formation and the ter- 
restrial planet formation stages that follow. Further out 
in the disk, this emission will control the photodissocia- 
tion rate of C2H and H 2 CO . which modify the molecu lar 
abundances at larger radii ([van Zadelhoff 'eTaL|[200l . 

Far-UV continuum emission has been considered be- 
fore in low-mass protoplanetary disk models, but some 
models assume that the A < 1400 A flux is dominated by 
the interstellar radiation field (ISRF) and stellar emis- 
sion lines. We present, however, evidence that the ab- 
solute flux of the far-UV continuum is several orders of 
magnitude greater than the ISFR. Integrating the con- 
tinuum fits (shown in Figure 3) over the far-UV bandpass 
(912 - 2000 A), we find that the absolute strength of the 
far-UV continua are 10 5-7 times the interstellar value 
(G = 1.6 x 10~ 3 erg cm~ 2 s -1 ; Habing 1968). Assum- 
ing distances of 140 and 70 pc for DF Tau and V4046 
Sgr, we find \og w {G/G ) = 6.8 and \og w {G/G ) = 5.6, 
respectively, at 1 AU. We note that this is only 912 - 
2000 A accretion continuum radiation field strength, not 
including the flux from atomic emission lines. Our radi- 
ation field estimates are consistent with the CTTS UV 
intensities calculated for 100 AU bv lBergin et al.l (|2004f ). 

One caveat is that this A-star-like accretion emission 
may not be present at all times, or at least not at the flux 
levels we observe. Figure 3 shows that optical activity 
may be a good proxy for the level of far-UV continuum. 
Thus the far-UV accretion emission is most likely vari- 
able, while the external field provided by the ISRF is 
constant. 



7 The fit to the DF Tau data reaches zero at 1122 A, but this 
is clearly an artificial cut-off imposed by the simple linear model. 
The data show appreciable continuum flux above the linear fit at 
A < 1200 A. 



5.2. Physical Conditions 0/ the Molecular Gas in the 
Inner Disk 

We find that the electron-impact excited H2 population 
in V4046 Sgr is characterized by N(R 2 ) ~ 10 18 cm" 2 
and T(H2) ~ 3000 K. This temperature is high but is 
not unreaso nable for the inner reg ion of an X-ray irra- 
diated disk HMeiierink et al.l 1200811 . These authors find 
that this temperature can be maintained in the inner 
« 20 AU of the disk, consistent with H2 em ission excited 
by ph otoelectrons produced by the X-rays (|Bergin et al.l 
2004). While the outer radius of the collisionally ex- 
cited H 2 may be as large as 20 AU, we suggest that this 
process could operate much closer to the central object, 
most likely < 1 AU (jlnglebv et a 1.1 12009:. The average 
width of the photo-excited H2 lines listed in Table 3 is 
52 ± 8 km s , with evidence of double peaked structure 
in the highest S/N lines. While some of this width may 
be due to turbulent motions within the disk, this sug- 
gests a concentration of photoexcited H2 at 0.19 AU, for 

1 — 35° and the simple assumption of a Keplerian disk 
where r H2 = (GM/< v FWH m >) 1 ^ 2 - 

If this is the inner edge of the warm molecular disk, it 
is co-spatial with the inner edge of the warm dust disk 
(0.18 AU; Jensen & Mathieu 1997), similar to the photo- 
excited H 2 distribution around the accreting brown dwarf 
2M1207 ([France et al.ll2010aft . While we cannot deter- 
mine the spatial origin of the electron-impact excited H2 
emission conclusively, we favor an origin in the ~ 0.2 - 

2 AU inner disk. While the electron-impact H2 emis- 
sion may be present at the disk midplane where plan- 
ets are actively forming, the observed H2 emission most 
likely originates between the disk midplane and the the 
disk surface. This radiation is being emitted from planet 
forming radii in the disk, but the collisionally excited 
H2 is not tracing the cold (T ~ 30 K), high column 
(N H > 10 22 cm -2 ), low electron fraction (x e < 10~ 6 ) 
disk midplane. 

The far-UV continuum and the "PDR-like" environ- 
ment of the inner disk call into question the assump- 
tion that non-thermal electrons are the only excitation 
source for the broad-band H2 emission that we observe 
in V4046 Sgr. We note that the best fit column den- 
sity and temperature values for the collisionally excited 
H2 are similar to those found for the Lya-pumped H2 in 
other CTTS disks [T(H 2 ) ~ 2500 K, iV(H 2 ) ~ 10 18 ~ 19 
cm~ 2 ] (Ardila et al. 2002; Herczeg et al. 2004), even 
though these states are not fully thermalized. At moder- 
ate spectra resolution, the spectrum of H2 excited by the 
far-UV continuum of a hot star is qualitatively very sim- 
ilar to the far-UV electron-impact spectrum (compare, 
e.g., Witt et al. 1989 and France et al. 2005 with Gustin 
et al. 2006). 

Finally, in §3.3.2 we suggested that hydrocarbons are 
responsible for the selective reddening of the observed H2 
electron-impact spectrum. Methane, the most important 
opacity source for electron-impact H2 emission in giant 
planet atmospheres, is difficult to observe at radio wave- 
lengths owing to its lack of permanent dipole moment. 
However, methane has been observed in higher mass 
prot ostellar disks at abu ndances w 10~ 6 that of hydro- 
gen (iBoogert et aLll2004D . Furthermore. iCarr k, Naiital 
(2008) presented infrared spectra of acetylene emission 
in the disk of the CTTS AA Tau. They posited that 
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the large observed abundance of C2H2 is a sign of ver- 
tical mixing within the disk, bringing the hydrocarbons 
up from the midplane. An analogous scenario fits our 
methane assumption. If we assume that the CH4 is 
~ 10~ 6 times as abundant as hydrogen in the planet 
forming region (Njj ~ 10 22 ~ 23 cm~ 2 ) and that this gas 
can be transported vertically towards the disk surface, 
then A^CH/t) is sufficient to attenuate the discrete H2 
line emission in V4046 Sgr [iV(CH 4 ) > 6 x 10 16 cm" 2 ]. 
Thus, we may have found indirect evidence for a reser- 
voir of hydrocarbons in the inner disk, where they are 
expected to become impor tant components of ex trasolar 
giant planet atmospheres (jSudarsky et al.l [20031 ). 



continuum can be described by a simple linear fit, which 
connects to the optical Balmer continuum given assump- 
tions about the optical activity level at the time of the 
far-UV observations. Accretion introduces an additional 
source of far-UV continuum emission that is 5 - 7 or- 
ders of magnitude stronger than the ISRF and has an 
effective temperature of ~ 10 4 K in DF Tau. We also 
find that the electron-impact H2 spectrum in V4046 Sgr 
is characterized by T(H 2 ) = 3000 K and N(H 2 ) ~ 10 18 
cm~ 2 , with E e ~ 50 - 100 eV. Finally, we propose that 
hydrocarbon absorption is the most likely explanation for 
the differential reddening inferred from the attenuation 
of collisionally produced H2 emission lines below 1400 A. 



6. SUMMARY 

We have presented here the highest-quality observa- 
tions of the far-UV continuum in CTTSs acquired to 
date. Because the gas-rich CTTS phase is roughly co- 
temporal with the epoch of giant planet formation, our 
observations represent an important portion of the ener- 
getic radiation environment present during the formation 
of exoplanetary systems. We characterize the spectra of 
hot accretion emission and electron-excited H2 in two 
prototypical CTTSs. We find that the far-UV accretion 
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